Bound states in the continuum in open Aharonov-Bohm rings 
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Using formalism of effective Hamiltonian we consider bound states in continuum (B1C). They are 
those eigen states of non-hermitian effective Hamiltonian which have real eigen values. It is shown 
that BICs are orthogonal to open channels of the leads, i.e. disconnected from the continuum. 
As a result BICs can be superposed to transport solution with arbitrary coefficient and exist in 
propagation band. The one- dimensional Aharonov-Bohm rings that are opened by attaching single- 
channel leads to them allow exact consideration of BICs. BICs occur at discrete values of energy 
and magnetic flux however it's realization strongly depend on a way to the BICs point. 
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I. INTRODUCTION. 

In 1929, von Neumann and Wigner [| firstly pointed 
to the existence of discrete solutions of the single-particle 
Schrodinger equation embedded in the continuum of pos- 
itive energy states. Their analysis has been examined 
by Stillinger and Herrick in the context of possible 
bound states (BICs) in atoms and molecules. It has 
been demonstrated by Newton Q that strong coupling 
between scattering channels can give rise to BIC. BIC 
can be observed in the stationary transmission as reso- 
nant states with width which tends to zero as at least two 
physical parameters vary continuously as it was formu- 
lated by Friedrich and Wintgen 4] , who have also given 
the example of the hydrogen atom in a magnetic field. 
Such a BIC is a very fragile structure. A small perturba- 
tion transforms it into narrow resonance. Nevertheless, 
Capasso et al Q have reported direct evidence for BICs 
in a semiconductor superlattice. 

For better understanding of the phenomenon of BICs 
in transport through electronic devices it is useful to 
study as simple quantum system as possible. Robnik 
has shown that a simple separable two-dimensional 
Hamiltonian can develop BIC under perturbation of open 
channels. An explicit proof of an existence of BICs was 
presented recently by Cederbaum et al Q in the molec- 
ular system, if the electronic and the nuclear motions 
are coupled. In the present letter we consider the open 
Aharonov-Bohm (AB) rings which are good candidates 
to observe BICs for the external magnetic field and en- 
ergy of incident electron can be easily varied experimen- 
tally. Moreover the one-dimensional AB rings allow to 
treat BICs wholly analytically. A phenomenon of zero 
resonance widths at discrete values of energy of incident 
particle and some relevant physical parameter was estab- 
lished in many works since the work by Shahbazyan and 
M.E. Raikh |E 0, E3, IH 13 EE E3, llE HI, amongof 
them external magnetic flux was considered in |l7l Il8| . 
In this letter we focus on the scattering wave function in 
the vicinity and at BICs point and how BIC participated 
in transport. 



II. THE ONE-DIMENSIONAL RING 

Following Xia [Tgj we write the wave functions in the 
segments of the structure shown in inset of Fig. ^as 



ipi(x) = exjp(ikx) + r exp(— ikx), 

i>2{%) = a\ exp(ifc~a;) + a 2 exp(— ik + x), 

i^six) — h\ exp(jfc + a;) + b 2 exp(— ik~x), 

ip4(x) = teKp(ikx), 



(1) 



where k~ — ft — 7, k + = k + 7, 7 = 27r$/$ , $ = 
BirR 2 is the magnetic flux, $0 = 2irhc/e. The ring length 
2irR is chosen as unit. The boundary conditions (the 
continuity of the wave functions and the conservation of 
the current density) allow to find all coefficients in QJ. 
We write the corresponding equation in matrix form 



F1> = ?, 
where F(k, 7) is the following matrix 



(2) 



/-I 

-1 






-1 

-1 

1 



1 



,ik~ /2 





V -1 ^ 



1 



e -ik+/2 


k 4 

fc+ -iii 
~TT e 2 




1 



e »fc+/2 
fe+ 

k 4- 

k + iii 

V e 2 




1 



D -ik~ /2 



(3) 

cf T = (1 1 1 0). The vector ^ = (r t 01 a 2 h b 2 ) 
is the solution for the scattering wave function: 



r = 2(3cosfc-4cos7 + l)/Z, 

t = 16i(sin|cos2)/^ 

01 = 2(2^ - 3e-^ + 1)/Z, 

a 2 = 2(e ik + 1 - 2e^)/Z, 

Z = 8cos 7 -9e~ ifc -e ifc + 2, 



(4) 



bi, 2(^,7) == a i. 2(^,-7)- In Fig. [I]we show lines of the 
transmission zeros (\t(k, 7)] = 0, dashed lines) which 
cross the lines of the transmission ones (\t(k, 7)! = 1, 
solid lines) at points 



k m — 2"Km, m = ±1, ±2, . 
7„ = 27rn, n = 0, ±1, ±2, . 



(5) 
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As can be seen from the expression for the denominator 
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FIG. 1: Transmission zeros \t\ 2 — and ones \t\' 2 = 1 of the 
one-dimensional ring as function of the wave number k and 
flux 7. The zeros (ones) are shown by dashed (solid) lines. 
The thin solid lines represent the eigenenergies of the closed 
ring. 

Z in Eqs. the imaginary part of the poles vanishes 
at these points. Simultaneously, there is a degeneracy 
of eigenenergies of closed ring (k m — j) 2 . Here m is the 
magnetic quantum number that defines the eigen func- 
tions of the closed ring ip m (x) = e~xjp(ik m x) . The point 
k = is excluded from the consideration since it gives 
zero conductance. The peculiar points JSJ were shown in 
(l7| for the case of single lead attached to the Id ring. To 
show that the bound states in continuum (BICs) appear 
at the points (J5J, let us consider one of the points, say, 
s o = (fcii 7i) = 27r(l, 1). All the other points are equiv- 
alent because of the periodical dependence of the system 
on k and 7. In the vicinity of the point Sq we write Eqs 
Q in the following approximated form 



t 



Afc 



Afc 



«i 



i(A 7 ) 2 /2' 
3Afc + 2A7 
4Afc + 2iAj 2 



0-2 



i(3Afc 2 - 4A 7 2 ) 
4(2Afc-HA 7 2 ) [ 
Afc - 2A7 
4Afc + 2iA 7 2 



(6) 



where Afc = fc — fci, A7 = 7 — 71. The transmission 
amplitude in the vicinity of the point So in © is similar 
to the expressions obtained for a shifted von Neumann- 
Wigner potential j^. One can see that all amplitudes 
<2i,2) ^1,2 °f t ne inner wave functions are singular at the 
point So. Such a result for the BIC points was firstly 
found by Pursey and Weber |2^| . 

Eqs (J2J and © allow to show that the point s corre- 
sponds to the BIC one in an open one-dimensional ring. 
At this point the matrix © takes the following form 
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The determinant of the matrix F(sq) equals zero. There- 
fore, Ff = 0. By direct substitution of the vector 

fo =|(0 01— 1— 11) one can verify that / is 
the right eigenvector which is the null vector. The corre- 
sponding left null eigenvector is / = ^(— 111 — 10 0). 
It is well known from linear algebra, that if the determi- 
nant of matrix F is equaled to zero, then the necessary 
and sufficient condition for existence of solution of the 
equation is that the vector is orthogonal to vec- 
tor if [24j. In holds, indeed, / • if = 0. The solution 
of Eq. J2J at the point Sq can therefore be presented as 



(8) 



where a is an arbitrary coefficient and ip p is particular 
transport solution of Eq. (J2J) . By direct substitution one 

can verify that ip p =(0l|||^)is the particular so- 
lution of Eq. J2J. It is worthwhile to note that this result 
completely agrees with the scattering theory on graphs 
pfl l2l| . Texier [2(j has shown that for certain graphs 
the stationary scattering state gives the solution of the 
Schrodinger equation for the continuum spectrum apart 
for discrete set of energies where some additional states 
are localized in the graph and thus are not probing by 
scattering, leading to the failure of the state counting 
method from the scattering. Similar result was indepen- 
dently obtained recently by Voo and Chu [2^ • 

In the vicinity of the BIC point so the scattering state 
using 10 becomes, to leading order of Afc, A7, 



V>(s) 



A 7 7 + Afc^ p 
Afc-HA 7 2 /2 ' 



(9) 



where s = (fc,7). Thus, the scattering state in the near- 
est vicinity of the BIC point also is superposed of the 
BIC vector ^ and of the particular solution tp p . Eq. 

© shows that the limiting scattering wave state ip de- 
pends on a way s — > So If we at first take A7 = 0, then 
obtain %j> = ip which is a transport solution. If we, how- 
ever, choose at first Afc = 0, then have = j^~T i-c. 
the scattering state is diverging interior the ring. This 
formula shows that the BIC state can ^ e extracted 
from the scattering state by a special limit in ©. 



III. TWO-DIMENSIONAL DEVICES 

Typical open two-dimensional structures are dots or 
rings with attached leads. Numerically the transmission 
through them are solving by finite-difference equations 
which are equivalent to the tight-binding lattice model 
l25l l2a . The case of the quantum dots was considered 
in [l6|. Here we present results of computation for the 
two-dimensional ring with symmetrically attached iden- 
tical leads. In Fig. [21 we show the transmission zeros 
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FIG. 2: Zeros (dashed blue lines) and ones (solid red lines) 
of the transmission probability of the two-dimensional ring 
as function of wave number k = y/E — n 2 and flux 7 = 
BttR 2 /$o, $0 = 2-Khc/e. R — 2.5 is mean radius of the 
ring. The width of the ring and those of the leads are equaled 
to unit. The eigenenergies of the closed two-dimensional ring 
are shown by thin lines. The BIC points are marked by open 
circles and star. 



(dashed lines) and the transmission ones (solid lines) for 
the single-channel transmission. 

In order to find the positions and widths of the 
resonance states, we explore the non hermitian effec- 
tive Hamiltonian, which can be obtained by projection 



of whol e sy stem onto Hilbert space of closed system 
[2rj l27l |2S| . The effective Hamiltonian in the basis of 
closed system's eigenvectors can be written as [2F 

(b\H eff \b')=E b 5 bb ,-J2 



E 

C=L,R 



V h c VS e ik ? 



26] 



(10) 



Here E b and \b) are the eigenvalues and the eigenfunc- 
tions of closed system are given by quantum numbers b, C 
enumerates left and right leads and p does the open chan- 
nels of leads. A recipe to calculate the matrix elements 
V b c is given in |2y, |27| . Because of energy dependence of 
the effective Hamiltonian the positions and widths of the 
resonance states arc defined by the following nonlinear 
fixed point equations 

E x =Re(zxh,E x )), 2T X = -/m(z A ( 7 , E x )). (11) 

Here z x are the complex eigenvalues of the effective 
Hamiltonian I|1U|) H e ff\X) = zx\X) with right eigenstates 
I A). All the points at which Fx = 0, i.e. the width of the 
resonant transmission vanishes are marked in Fig. [21 by 
open circles and star. 

The equation for the scattering wave function mapped 
interior the ring \ipR) can be derived from the Lippmann- 
Schwinger equation [26l I2H |28| and takes the following 
form 



Here V L is the coupling matrix between the left lead 
and the ring provided that a particle incidents from the 
left lead in the first channel. This formula is similar to 
© for the Id ring. If Det(H eff - E) ^ 0, then in the 
biorthogonal basis | A) the scattering wave function takes 
a simple form pfl |28| 



VP 1 



E 



E-z x {7,E) 



|A), 



(13) 



where 

Vx = {X\V\E,L,p=l) = J dy B ^x(yB)sm(7ry B ), (14) 

tpx are the left eigen functions of H e ff, yg runs over the 
boundary that connects the closed ring and the left lead 
with the first channel excited {p — 1). We assume that 
magnetic field subjects only the ring. 

Let us denote a set of physical parameters of the sys- 
tem as s. For example, for present case of the ring 
s = (_E,7), although for the quantum dot s might be 
energy and confined potential [l|j. Let us consider 
the point So = (i?o,7o) at which Eq. Ijllfl i s fulfilled 
Eo = za (-Eo,7o) and Tx„ = 0, i.e. one of the complex 
eigenvalues of H e ff is real at this point. For E = Eq one 
have equality (H e fj — E)\X ) = 0. Comparing this equa- 
tion to l|12|l we see that the eigen state |Ao) corresponds 
to the solution of the Lippmann-Schwinger equation if 
there were no ingoing current in the left lead. Respec- 
tively, the state |Ao) can not give rise to outgoing currents 
because of the continuity equation for the current den- 
sity. In order to fulfill that we have to consider that the 
eigen function ipx does not overlap with the first channel 
of the left lead, i. e. V^(so) = 0. This may be also es- 
tablished by consideration of the transmission amplitude 



(E,L\V L \X)(X\V R \E,R) 
E-zx 



(15) 



Because of symmetry of the system relative to the left 



and right leads 



1^1 



In approaching the point 



s — > So the denominator E — Zx (s) — > 0. Therefore, in 
order the ratio \V x L {s)\ 2 / (E — zx {s)) remained finite in 
H15JI it is necessary |V^(s)| — > for s — > So. Thus, at the 
BIC point we have orthogonality of the righthand state 
(V\E,L,p = 1)) in Eq. JT3J| to the left eigen state (A |. 
Then, in full correspondence to the consideration of the 
Id ring (Eq. JHJ), we have the following solution for the 
scattering state interior the ring 



l^fl(so)) = a|Ao(so)) + \ip p (s )), 



(16) 



(H eff -E)\^ 



V L \E,L,p=l) 



(12) 



where coefficient a is arbitrary. Right eigen function 
ipx (s ) °f the effective Hamiltonian is squared integrable 
and therefore is the BIC function shown in Fig. [3] (a). 
Although the BIC function is disconnected from the the 
first channel of the left lead, it couples with the next 
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FIG. 3: The BIC function \lp\ 1 which is the eigen function of 
the effective Hamiltonian 1101 (a) and the transport solution 
\%p p \ (b) at BIC point marked in Fig. |21by star. 

channels p > 1 of the leads which are evanescent modes. 
As a result the BIC function has exponentially small tails 
in the leads as might be seen from Fig. [3] (a). 

Moreover the coupling of the 2d ring with the evanes- 
cent modes gives rise to that as Fig. shows the BIC 
points are close to but different from points at which 
two eigen functions of closed 2d ring classified by mag- 
netic quantum numbers m have the same energy. The 
evanescent modes have imaginary wave numbers k p which 
change effectively the Hamiltonian of closed ring by ma- 
trix 

E E V b c p V b ^ p e-^~{d/RY 

p^l G=L,R 

via Eq. (|1(J|> . Therefore only for limiting case of the Id 
ring d/R — > the BIC state will consists of a pair of 
eigen states of closed ring as seen from Fig. ^and as was 
confirmed by computations. 

In the vicinity of So a value E — z\ (E,~/) is small. 
Then we can split the summation over A in (|f 3|) by two 
parts, A = Ao and A ^ Xq and similar to © write the 
scattering state as 

|^(s)) = a(s)|Ao(s)) + |V P (s)) ) (17) 

where 



and \ip p ) is contribution of all other resonances. 

As different from the Id ring, we can study behavior of 
singular coefficient (|18|) only numerically. Let us encircle 
the BIC point k — ^JE — n 2 , 70 as Afc — r cos 0, A7 = 
r sin (j) as shown in Fig. 0] (a) where r is the radius of 
encircling. Angular behaviors of quantities defining the 
parameter a are shown in Fig. 0](b,c). In particular the 
numeric shows that |Vb(s)| ~ |s — Sol 1 / 2 . The behavior of 
a in Fig. 01 (d) is very similar to the behavior of the pa- 
rameter a = A7/ (Afc + iA7 2 /2) for the Id ring (see ©) 
except that in the 2d ring we observe phase difference. 
As one can see from Fig. 0] (d) at <j) = 0, 7r ( A7 = 0) the 
parameter a — 0, and at 4> = (j)Q < 7r/2, a — > 00. The 
angle 4>q exactly corresponds to the direction of tangen- 
tial line of transmission zero shown by dashed line in Fig. 




4>/2tc §/2x 



FIG. 4: Angular behavior of quantities defining the parameter 
l|I8[l around the BIC point marked by star in Fig. 

121 Therefore in order to extract the \ip p ) from the scat- 
tering wave function 1)13(1 we should put at first A7 = 
and then limit Afc — > 0. If take limit to the BIC point 
along A7 — tan((f>o)Ak, the scattering state transforms 
to the BIC state | Ao) shown in Fig. 13(a). The particular 
solution for the scattering wave function \ip„) is shown in 
Fig. E|(b). 



IV. CONCLUSION 

Formulas (JSJ) and l|17|) are the key ones which show that 
scattering wave function ip is not unique since BIC can 
be superposed with arbitrary coefficient a to ip . Such 
kind of decomposition was established recently for the 
scattering theory on graphs [20L l2lj. Thus, at the point 
the So, the system becomes degenerate. Usual trans- 
port solution with energy E = Eq is complemented by 
the squared integrable (localized interior the ring) state 
|Ao(so)) with the same energy Eq orthogonal to the for- 
mer. The last state is therefore BIC. Our consideration 
shows exactly that BIC is the eigen vector of the non- 
hermitian effective Hamiltonian H e ff at those point at 
which the complex eigenvalue of H e ff becomes real and 
coincides with the energy of incident particle. The scat- 
tering matrix is unique but not analytical at the BIC 
points as could be seen from formula JSJ. As seen from 
there, the transmission zeros cross the transmission ones 
at the BIC point. Note that these results are not re- 
stricted by only AB rings but applicable for any open 
quantum system which allow to vary at least two physi- 
cal relevant parameters, for example, energy of incident 
particles and shape of billiard ^(| . 

BIC is disconnected from both single channel continua. 
In order to achieve that BIC is to be a such superposi- 
tion of eigen states of closed system that overlapping l(T4"|) 
vanishes at BIC points. Specifically in the present case 
of the AB ring attached to the single channel leads this 
superposition becomes odd function relative to the even 
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function of the leads as seen from Fig. ?? (a). For the Id 
ring nodes of BIC are to be at points of connection of the 
ring to leads, thereby at those poin ts where the ratio of 
lengths of the arms is rational [2(j. However for the 2d 
ring the leads are to be attached exactly symmetrically 
as shown in Fig. ??. It follows then that a violation of 
symmetry of the system relative to transport axis x leads 
to breakdown of BIC. In particular it occurs for system 
disordered by impurities. In order BIC could survive 
under this violation of symmetry one can use geometry 
given in |29| in which infinite strip attached to the ring . 
Moreover impurities lift a degeneracy of closed ring [3(| . 
However as shown in 0] a condition for BIC to survive 
is still remaining. From above it follows that for the sys- 
tem symmetrical relative the y-axis (axis perpendicular 
to the transport axis) all odd eigen states of closed sys- 
tem are BIC's provided that the leads are excited in the 
first even channel. Then a perturbation which lifts this 
symmetry transforms BIC's into resonance states widths 
of which are proportional to the perturbation. The ex- 
ternal magnetic field which subjects only the ring is an 
example of such a perturbation. 

The electron-electron interactions preserve degeneracy 
of closed ring |30| . They modify the energy spectrum and 



the coupling between leads and closed ring. As shown 
in |13| variation of the coupling changes a position of 
BIC So however is not important to achieve real value 
of the complex eigen value of the effective Hamiltonian. 
However the Coulomb interactions might be important in 
respect that BIC's can exhibit discrete charging similar 
to that predicted for resonance trapping in quantum dots 
strongly connected to the leads |3l| . The strong coupling 
of closed quantum system with leads (\Vf p \ ^> \Eb — Ey \ 
in terms of fXD)) ) is hardly achievable while an existence 
of BICs is free of a value of the coupling between the 
closed system and continua. 

Processes of inelastic scattering give rise to finite res- 
onance width. In that sense BIC is very subtle phe- 
nomenon for electron transmission. However as formulas 
© and Ijl7(l show, BIC state participates in the scat- 
tering wave function. If above mentioned processes are 
efficiently small BIC state can dominate in the vicinity of 
the BIC point for proper choice of physical parameters, 
energy and flux, as shown in Fig. 01 (d) . 
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